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In this paper we study the two-phase size-structured population model with infinite states-at-birth and 
distributed delay in birth process. The model distinguishes individuals by two different status: the 'repro- 
ductive' stage and the 'nonreproductive' stage. We establish the well-posedness for this model and show 
that the solution of this model exhibits asynchronous exponential growth by means of semigroups. We 
also consider a special case in which the individuals in the 'reproductive' stage and the 'nonreproductive' 
stage have the same growth rates and give a comparison between this two-phase model with the classical 
one-phase model. 
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1. Introduction 

Among those size- structured population models with infinite states-at-birth, the fundamental one 
is as follows: 

dn d(y(x)n) f~ a 

H = -fi(x)n(t 9 x) + / P(x,y)n(t,y)dy, 0 < x < a, t > 0, 



dt dx 



n(t,0)=0, t>0, (1) 
n(0,x) = hq(x), 0 < x < a. 

Here, the unknown function n(t,x) denotes the density of individuals of size x e [0,a] at time 
t e [0, oo), where a > 0 represents the (finite) maximum size of any individual in the population. 
The vital rates fi(x) and y(x) denote the death and growth rates, respectively. It is assumed 
that individuals may have different sizes at birth and therefore ft (x, y) denotes the probability 
which individuals of size y give birth to the individuals of size x. The asymptotic behaviour of 
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its solution can be obtained by using the generalized relative entropy method [12]. Recently the 
similar nonlinear model was studied in [10,11]. An interesting situation to consider is when the 
birth process undergoes an observable period so that there is a time lag between conception and 
birth [6,15,16,18]. In such considerations, the model (1) should be improved by the following 
model with the distributed delay: 

dn d(y(x)n) f~ a f° 

1 = -fi(x)n(t,x) + / / P(a,x,y)n(t + a,v)dady, 0 < x < a, t > 0, 

dt dx J 0 J_ T 

n(t, 0) = 0, t > 0, 

n(cr,x) = h(cr,x), a e [— t,0], 0 < x < a. 

(2) 

Here (a , x, y) denotes the rate at which individuals of size y give birth to the individuals of size x 
after a time lag —a starting from conception and r is a constant denoting the maximal delay. For 
example, in a host-parasite system, the distributed delay in model (2) may be given by the time 
lag between laying and hatching of the parasite eggs [2]. Moreover, unlike the non-distributed 
delay case, the time lag considered here can change from 0 to r , i.e. it is distributed in the interval 
[0, r]. The idea of considering the distributed delay is inspired by the work of Piazzer and Tonetto 
[14], where a different age- structured population model with distributed delayed birth process 
was studied. Recently, the model (2) has been proved in [4] that the problem is globally well- 
posed, and the solutions exhibit so-called asynchronous exponential growth (we refer the readers 
to [1,4,9,14] for definition). 

In this paper we study a model which distinguishes individuals by two different status: the 
'reproductive' stage and the 'nonreproductive' stage. In this model, only individuals in the 'repro- 
ductive' stage reproduce. In fact the 'reproductive' stage and the 'nonreproductive' stage occur 
in the evolution of many populations generally. We denote by p(t,x) and n(t,x) the densities of 
individuals in the 'reproductive' stage and the 'nonreproductive' stage of size x e [0, a] at time 
t g [0, oo), respectively. Then the model reads as follows: 

dp d(y\(x)p) 

— H = -ViWP ~ Pi(x)p + fawn 



dt dx 



pa pO 

+ vj J P(cr,x,y)p(t + cr,y)dcrdy for0<x<a, t>0, 



dn d(y 2 (x)n) 

H « = —H>2\x)n + PiWP - Pi(x)n 



dt dx 



(3) 

+ (1 - v) / / P(cr,x,y)p(t + a,v)dadv for0<x<a, t>0, 



pa pO 



p(f,0) = 0 forr>0, 
n(t, 0) = 0 for t > 0, 

p(cr,x) = p(cr,x) for a G [— r,0], 0 < x < a, 
n(0,x) = rio(x) for 0 < x < a. 

Here y\(x) and 72 CO represent the growth rates of the individuals in the 'reproductive' stage 
and in the 'nonreproductive' stage, pi(x) and P2CO represent the transferring rates between the 
'reproductive' stage and the 'nonreproductive' stage, respectively, and fii(x) and /X2CO represent 
the death rates of the individuals in the 'reproductive' stage and the 'nonreproductive' stage, 
respectively. Also fi(cr,x, y) represent the rate at which the individuals in the 'reproductive' stage 
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of size y give birth to the individuals in the 'reproductive' stage or the 'non-reproductive' stage of 
size x after a time lag —a starting from conception, and v is a constant, 0 < v < 1. In addition, 
p and no are given functions defined in [— r, 0] x [0, a] and [0, a], respectively. Later on we shall 
denote 

p Q (x) = p(0,x) for 0 < x < a. (4) 

The similar model without distributed delay has been proved in [9] that the problem is globally 
well-posed, and the solutions exhibit asynchronous exponential growth. The purpose of this work 
is to extend the results in [4,9] to the model. We shall prove that under suitable assumptions on 
fi\(x), /X2(X), Pi W, Pi(x), Yi(x), and (p, no) the model (3) is globally well-posed, 

and its solution possesses the properties of asynchronous exponential growth. 

Throughout this paper, ^(x), P\(x), P2 (x), y\(x), Y2(x) and ft(cr,x, y) are supposed to 

satisfy the following conditions: 

(H.l) fix, /JL2, pi and p 2 are nonnegative and continuous functions defined on [0, a]. 
(H.2) y u y 2 e C l [0,a] and y x (x), y 2 (x) > 0 for all* e [0,a]. 
(H.3) 0 e C([-r,0] x [0,0] x [0,5]) and£ > 0. 

In order to prove the property of asynchronous exponential growth, we make the additional 
assumptions: 

(H.4) If 0 < v < 1, p 2 (» > 0 for all x e [0, a]. If v = 1, pi (jc), p 2 (x) > 0 for all x e [0, a]. 
(H.5) y) > Ofory > *. 

We introduce the subspace ^4 of W l,l ((—T,0),L l [Q,a\) and the subspace S of W l,l (0,a) as 
follows: 

^4 : = {p e W hl ((-r,0),L l [0,a\) : p 0 e W hl (0,a),p 0 (0) = 0}, 
B := {n 0 e W 1J (0,a) : n 0 (0) = 0}, 

where £ 0 is defined by Equation (4). Note that since W hl ((-r , 0), L 1 [0, 5]) c C([-r , 0], L 1 [0, a]), 
for any £ e W l > l ((-r,0),L l [0, a] ),p 0 is well denned. Similarly, since W l > l (0,a) C C[0,a] ?j po(0) 
and ^o(0) are also well defined. 

Our first main result establishes the global well-posedness of the problem (1) and reads as 
follows: 

Theorem 1.1 For any (p,no)tAxB, the model (3) has a unique solution (p,n) e 
(C([-r, 00), L 1 [0, a]) n C([0, 00), W 1 * 1 (0, a)) n C 1 ([0, 00), L 1 [0, a])) x (C([0, 00), 
W u (0,a)) n (^([O, 00), L 1 ^),^])). Moreover, for any T > 0, f/ie mapping (p,n 0 ) i-* 
(p,n) from A x B to (C([-r, oo),/, 1 ^, 5]) n C([0, 00), W u (0, a)) n ^([0, oo),/, 1 ^, a])) x 
(C([0, 00), W u (0,a)) n C^tO, 00),^ [0,5])) w continuous. 

The proof of this result will be given in Section 2. 

Actually, from the proof of Theorem 1.1 we shall see that for any F e £ , where 

£ := L l ([-r,0],L l [0,a\) x L l [0,a], 

the model (3) has a unique so-called mild solution (p, n) in the sense that [t i-> [(cr,x) i-> 
p(t + cr,x)]] e C([0,oo),L l ([-r,0],L l [0,a]), [t \-> [x i-> n(t,x)]] e C([0, 00), L 1 ®, a]) and it 
satisfies an integral equation which is equivalent to Equation (3) in a suitable sense. This defines, 
for each t > 0, an operator Tit) : £ —> £ , i.e. T(t)F := + •, -),^(^ •))• Later, we shall see 
that (T(0)*>o is a strongly continuous semigroup in £ . 
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Our second main result of this article studies the asymptotic behaviour of this semigroup and 
reads as follows: 

Theorem 1.2 There exist a rank one projection n on £ and constants Xq g R, s > 0, M > 1 
such that 

\\Q- Xot T(t) -Tl\\<M Q~ £t for all t > 0, (5) 
where || • || denotes the operator norm on £. 

The proof of this result is given in Section 3. The parameter Ao is called intrinsic rate of natural 
increase or Malthusian parameter [17]. Theorem 1.2 shows that the solutions of the model (3) 
exhibit asynchronous exponential growth. 

Next we consider the special case of the model (3), where y\ (x) = y 2 (x) = y (x) and give a com- 
parison between this two-phase model with the one-phase model. We want to give a comparison 
between the asymptotic behaviours of the sum of the densities of individuals in the 'reproductive' 
stage and the 'nonreproductive' stage and the solution of the one-phase model after modifying 
the death rates and fi(cr,c,y) properly. Note that the above result says that there exists a positive 
vector function (w, v) g £, such that for any (p, no) G A x B, the mild solution (p(t + •,•)> n(t, •)) 
of the model has the following asymptotic expression: 

(p(t + ■,■), n(t, •)) = [c(w(-, •), v(-)) + 0(e~ £t )] e^ as oo, 

where c is a constant uniquely determined by the initial data (p, no); c > 0 provided (p, no) > 0 
(i.e. (p,no) > (0,0) and (p,n 0 ) # (0, 0)). We denote 

N(t,x) = p(t,x) + n(t,x) for t > 0, 
N(cr,x) = p(cr,x) + no(x) for cr G [— t, 0]. 

We want to compare N(t,x) with the solution N(t,x) of the problem 

dN d(y(x)N) - r f° 

9* 9* Jo J-t 

^,0) = 0, t>0, 

N(cr,x) = p(cr,x) +no(x), a G [— r,0], 0 < x < a, (6) 

where 

w(0,jc) 

C/(jc) = w(0,x) + v(x), 6»(x) = 2 \ 

U(x) 

fi(x) = 0(x)tn(x) + (1 - 0(x))fjL 2 (x), Pi(cr,x,y) = 0(x)P(cr,x,y). 

We can see that 6(x) is the asymptotic proportion of the individuals in the 'reproductive' stage 
in the population. Since the model (6) describes the evolution of the sum of the densities of 
individuals in the 'reproductive' stage and the 'nonreproductive' stage in the asymptotic sense, 
one might expect that N(t,x) — N(t,x) —> 0 as t —> oo. But to our surprise, this is actually not 
the case. In fact, we have the following result: 

Theorem 1.3 Let the notation be as above. We have the following relation: 

N(t, x) Q~ Xot - N(t, x) Q- Xot = cll(x) + 0(Q~ £t ) as t -> oo, 
where c is a constant which is generally non-vanishing. 
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The proof of this result will be given in Section 3. Theorem 1.3 shows that the asymptotic 
behaviours of the sum of the densities of individuals in the 'reproductive' stage and the 'nonre- 
productive' stage and the solution of the one-phase model are different and the research of the 
model with two stages is meaningful. 

The layout of the rest of the paper is as follows. In Section 2 we reduce the model (3) into an 
abstract Cauchy problem and establish the well-posedness of it by means of strongly continuous 
semigroups. In Section 3, we prove that the solution of the model (3) has asynchronous exponential 
growth. In Section 4, we give the proof of Theorem 1.3. 



2. Reduction and well-posedness 

In this section we reduce the problem (1) into an abstract Cauchy problem and establish the well- 
posedness of it by means of strongly continuous semigroups. We refer the reader to see [5,14] for 
similar reductions. 

First, we introduce the following operators on the Banach spaces X := L l [0,a] x L l [0,a], 
Y 0 := {u e W hl (0,a) : u(0) = 0} and E := L l ([-r, 0],L l [0, a]): 

A( W ,v) = (-(n(-)^) / ,-(K2(-)v) / ), 
with domain D(A) = Yo x Yq, 

B(u, v) = (-0*1 (•) + Pi(-))u + p 2 (-)v, Pi(-)u - (p 2 (-) + M2(0)v) for (u, v) e X, 
C(u) = (vCi(u),(l — v)C\(u)) for u e E, 

where 



pa pu 

C 1 (u) = J j P(cf,>,y)u(t + cr,y)dady. 



We note that A e C(D(A),X), B e C(X) and C e C(E,X). Using this notation, we rewrite the 
model (1) into the following abstract initial value problem for a retarded differential equation in 
the Banach space X\ 

d(u(t),v(t)) 

' = (A + B)(u(t), v(0) + Cu t , t > 0, 
at 

(w(0),v(0)) = (po,no), (7) 

U 0 = p, 

where u : [0, +oo) —> L l [0, a] and v : [0, +oo) —> L l [0, a] are defined as u(t) := p(t, •) and 
v(t) := n(t,-), respectively, and u t : [— r,0] —> L l [0,a] is defined as u t (cr) := p(t + cr),cr e 
[-r,0]. 

Remark 2.1 As usual, if the functions u : [— r, oo) -> L l [0, a] and v : [0, oo) -> L l [0,a] sat- 
isfy u e C([-r,oc), W u (0,a)) nC([0,oc),y 0 ) nC l ([0,oo),L l [0,a])mdv e C([0,oo),7 0 ) n 
C 1 ([0, oo), L 1 [0, a]), respectively, and the problem (7) in usual sense, we say that functions (u, v) 
is a classical solution of the problem (7). It is evident that a necessary condition for the prob- 
lem (7) to have a classical solution is that p e W 1,1 ([—r,0],X) and the functions po defined by 
Equation (4) and no belong to Yq. 
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Next, we introduce the following operators in the Banach space E: 

(Gu)(a) := —u with domain D(G) = W u ([-r, 0],L l [0, a]), 
da 

Qu := u(0) foru e D(G). 

We note that G e £(D(G),E) and Q e £(D(G),X). We now let E:=£xl, and introduce 
operators Ao, B and A in E as follows: 



A °((w,v)) (o A + b)\(u,v) 



with domain D(A 0 ) := { ( f „ e D(G) x D(A) : Qu = u^, 

»(£>)=-(? o) («;:,„) H->) eE - 

A = A 0 + B with domain D(A) := D(A 0 ). 

We note that A 0 G £(D(A),E), B e £(E,E) and A e £(D(A),E). Using this notation, we see 
that the problem (7) can be equivalently rewritten into the following abstract initial value problem 
of an ordinary differential equation in the Banach space E: 

U'(r) = AU(0, t > 0, 

(8) 

U(0) = Uo, 



where 



U(,) =U(0?v(0)J ^ U °= («(0) M v(0)) 



Remark 2.2 As usual, we say that a function U : R+ := [0, oo) —> E is a classical solution of 
the problem (8) if U e C([0, oo),D(A)) n C^tO, oo),E) and if it satisfies Equation (8) in the 
usual sense. 

To be rigorous, we write down the following preliminary result: 

Lemma 2.1 Let the necessary condition mentioned in Remark 2.2 is satisfied. If (u(t),v(t)) 
(u : [— r, oo) -> L l [0, a], v : [0, oo) -> L^O, a]) w a classical solution of the problem (7), 
U(0 = ( ( u (t)'v(t)) ) ^ a classical solution of the problem (8). Conversely, if\J is a classical solution 
of the problem (8), then U has the form U(7) = ( (u(t)Xt)) ) / or a ^ t > 0> a/id Z?y extending the first 
component of its second component u = w(7) to [— r, oo) swc/z that u(t) = uofor t e [—z,0),we 
have that (u, v) is a classical solution of the problem (7). 

Proof See Lemma 2. 1 of Bai and Xu [4] and Theorem 2.2 of Piazzer and Tonetto [14]. ■ 

In the sequel, we consider the semigroup generated by the operator A. We first consider the one 
generated by the principle part Ao of A. We have the following results: 

Lemma 2.2 The component operator A + B of operator matrix Ao generates a strongly 
continuous semigroup {7b(f)}f>o on ^' 
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Proof We note that A generates a nilpotent semigroup on X(see Theorem 2.2 of Farkas and Hinow 
[9]). Since B e C(X), by using the perturbation theorem for generators of strongly continuous 
semigroups in Banach spaces (see Theorem III. 1.3 of Engel and Nagel [8]), we get this lemma. ■ 

Lemma 2.3 The operator A 0 generates a strongly continuous semigroup (T 0 (0)/>o on E, given 

where (S(t)) t >o is a nilpotent left shift semigroup on E, given by 
(S(t)f)(a,x) = 



f(t + cr,x) ifcr + t<0, 
0 ifcr + t>0. 



and T t : X —> E are linear operators defined as 



(10) 



\jri(T 0 (f + <T)F) ifa+t>0, 
T <^ F -=\ 0 ifa+t<0, (U) 



where it\ is the projection onto the first coordinate. 



Proof See [5, Theorem 2.2] and Lemma IV. 1 .2 of Engel and Nagel [8] . ■ 

Since B e £(E), by using the perturbation theorem for generators of strongly continuous 
semigroups in Banach spaces (see [8, Theorem III. 1.3]), we get the following lemma: 

Lemma 2.4 The operator A generates a strongly continuous semigroup (T(0)f>o on E. 

By using the theory of strongly continuous semigroups in Banach spaces, we get the following 
result: 

Theorem 2.5 For any given initial data Uo = ( ( M (o) M ,°v(0)) ) G D(A), the problem (8) has a unique 
solution U e C([0, +oo);D(A)) n C l ([0, +oo);E), given by 

U(t) = T(t) ((u(0)M0))) f ° rt ^°- (12) 



By Lemma 2.1 and Theorem 2.5, we see that Theorem 1.1 follows. 



3. Asynchronous exponential growth 

In this section we study the asymptotic behaviour of the solution of the problem (1). We shall 
prove that the semigroup (T(0)?>o has the property of asynchronous exponential growth on E. 
We denote by &> e ssW the essential growth bound of the semigroup (T(t)) t >o with generator A, 
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i.e. 




(13) 



cdq(A) the growth bound, i.e. 



co 0 (A)= lim ^ log H7X0II 



(14) 



and s(L) the spectral bound , i.e. 



s(A) = sup{ReA, : X e cr(A)}. 



If we prove that the semigroup (T(t)) t >o is an irreducible positive strongly continuous semigroup 
(we refer the readers to [7,8] for definition) satisfying the inequality &> e ss(A) < cdq(A), then by 
[7, Theorem 9.10 and Theorem 9.11], the semigroup (T(0)?>o has the property of asynchronous 
exponential growth on E. Thus, in the sequel we step-by-step prove the above assertions about 
the semigroup (T(t)) t >Q. 

Lemma 3.1 The semigroup (T(t)) t >o generated by A is positive and eventually compact (we 
refer the readers to [7,8] for definition). 

Proof Since B is a positive bounded linear operator in E, the positivity of (T(0)f>o follows if we 
prove that the semigroup (T 0 (0)?>o generated by A 0 is positive (see Corollary VI. 1.11 of Engel 
and Nagel [8]). Since (S(t)) t >o is positive, by the expression (9) of (To(0)?>o> we only need to 
prove that the semigroup (7q(0)?>o generated by A + B is positive. The positivity of A + B can be 
easily obtained by the characteristic method (see [3, Lemma 2.1] and [9, Theorem 2.2]). Therefore 
the positivity of (T(t)) t >o follows. From the proof of [9, Lemma 3.6], we have that 7b (0 = 0 for 
t > T, where V = max{/ 0 a d§/yi(§), f Q a ^y}- From Equation (11), we also see that T t (t) = 0 
for t > r + r. Hence, from Equations (9) and (10), we have that To(0 = 0 for t > T + r. This 
particularly implies that (To(0)f>o is compact for t > V + r. Thus, by Proposition III. 1.14 of 
Engel and Nagel [8], the eventual compactness of (T(t)) t >o follows if we prove that B is compact. 
We note that the only nonzero component operator of operator matrix B is C : E — >► X. We use 
the method which is similar to Lemma 3.6 in [9] to prove that C is compact. Hence the desired 
assertion follows. ■ 

By the eventual compactness of the semigroup (T(0)?>o an d Corollary IV.3.12 of Engel and 
Nagel [8], the following result holds. 

Lemma 3.2 s(A) = 000(A). 

In the proof of Lemma 3.1, we have that the semigroup T 0 (0 = 0 for t > V + r . Then by the 
definition (14), we have that &> e ss(Ao) = —00. Since B is compact on E, by Proposition 2.12 of 
Clement et al. [7], we have the following result: 

Lemma 3.3 co ess (A) = — 00. 

Lemma 3.4 The semigroup (T(0)?>o generated by A is irreducible. 

Proof By Lemma VI. 1.9 of Engel and Nagel [8] and Lemma 3.4 of Bai and Cui [3], if we 
prove that {R(X,A)¥, ^) > 0 for some A > 0 and all F e E and V ^ e E* such that F > 0 and 
^ > 0, then the desired assertion follows. Let it\ and 7x2 be the projections onto the first and 
the second coordinates, respectively. We will prove that 7t\(R(X, A)F)(cr,x) > 0 for almost all 
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(a, x) e [-r,0] x [0,5], 7r 1 (7T 2 (^(A,A)F))(x) > 0 and 7r 2 (7r 2 (^(A,A)F))(x) > 0 for almost all 
x e [0, a]. To this end, we first deduce an useful expression of R (A, A). For F = (f(a,x), F(x)) e 
E, let U = R(k,A)F. Then U satisfies the equation 

(A/-A)U = F. (15) 

By writing U = (u(a,x), U(x)) and F = (f(cr,x),F(x))(U = (u, v), F = (f,g)), we see that the 
above equation can be rewritten as follows: 

d 

Xu(a,x) u(a,x) =f(cr,x), — t < a < 0, 0 < x < a, 

da 

d 

Xu(x) + — {y\{x)u{x)) = f{x) - ijl\{x)u{x) - p\(x)u(x) + p 2 (x)v(x) 
dx 

+ v / / P(a,x, y)u(a,y) da dy, 0 < x < a, 



no 
P(a,x,y)u(a,y)dady, 0 



d (16) 
Xv(x) + — (y 2 (x)v(x)) = g(x) - fi 2 (x)v(x) + p\(x)u(x) - p 2 (x)v(x) 
dx 



no 
P(cr,x,y)u(cr,y) da dy, 0 < x < a, 



zi(0,x) = w(x), 0 < x < a, 
M (0) = 0, v(0) = 0. 



Then 



S(<r,*) = e ka u(x) + e ACT / t~^f{^x)d^ 



f 

J a 



, , r , , ^/;/_ U T ^(a,,, } ;)e^(^) + /;e-^,v)d§)dad } ; 

Jo £ U (s)yi(s) 
V(X) = (1 - V)£ 2 A^) / =r~7\ 7T ^ 



+ ^2xW / ^"T7 — 

Jo 



where £ u (jc) = e~ fo(( x +m( s )+Pi( s )+y{( s ))/n( s )) ds andE^C*) = e^^^^+^C^+^C^+y^^))/^^))^^ 
For each A g C, we define two operators M k and Nx on E as follows: 

f(cr,x) \ ( e^/iW \ n7 , 
(fi(x),f2(x))J \(hi(x 9 k) 9 h 2 (x 9 k))J 9 K } 

(fi(x)Mx))J V (gi(x 9 k),g 2 (x 9 k)) J 9 K } 
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where 



t , r , , , vfZf° T P(a,s,y)^f 1 (y)dady + p 2 (s)Ms) 

h\(x,X) = E lk (x) / — — - ds, 

E n (s)y l (s) 



f 

Jo 

f 

Jo 

gl (x 9 X)=E lk (x) f — «> i_ _ ds , 

Jo 



^ r , , , (1 - v) /; J_ r ^(a, e^/i W da + (.) 

h 2 (x,k) = E 2k (x) / — — - ds 

E 2k (s)y 2 (s) 

V fo f-r P(°>S> y) ^ fa y) ^ dy + /i (j) 



, ^ _ , , f x d - ^ g f-r P(°> fOO e AfT g e-*/(g,y) dg da dy +/ 2 W . 
Jo E 2k (s)y 2 (s) 

Since 

\\M k (f(a 9 x) 9 (f(x) 9 g(x))\\ E ^0 asA^+oc, (19) 

there exists A* > 0 such that \\M k || < 1 for A > A*. This implies that the inverse (I — M k )~ l exists 
and is a bounded operator for A > A*. From Equations (17)— (19), we see that the resolvent of A 
is given by 

oo 

R(k,A) = (I - M k )~ l N k = J2 M x N ^ for x > ( 2 °) 

n=0 

We first consider the case in which 0 < v < 1. Let 0 < (f(cr 9 x) 9 (f(x) 9 g(x))) G E and 
(f(a,x), (f(x), g(x))) / 0. Without loss of generality, we assume that/(a,x) > 0 for almost all 
(cr,x) G [ao,ai] x [*o,xi] with some — r < oo < &\ < 0 and 0 < xo < x\ < a. Then for almost 
allx G [0,a], 

7Ti[7r 2 (^(f, (f, £)))](*) > 0, 7t 2 [n 2 (N k (f, (f,g)))](x) > 0. 

This is for almost all (cr,x) e [— r, 0] x [0, a], 

*i(M k N k (f 9 (f 9 g)))(a 9 x) >0. 

From Equation (20), we obtain that 7Ti(7?(A, A)F)(a,x) > 0 for almost all (a, x) G [— r,0] x [0, a] 
and n 2 (R(k, A)F)(jc) > 0 for almost all x e [0,a]. If we assume that/(x) > 0 for almost all 
x G [xo,xi] with some 0 < xo < x\ < a. Then for almost all x G [xo, a\, 

7t 1 [7t 2 (N x (f 9 (f 9 g)))](x)>0. 

This further implies that for almost all x G [0, a], 

7ri[7r 2 (M k N k (f, (f,g)))](x) > 0, 7t 2 [7t 2 (M k N k (f 9 (f 9 g)))](x) > 0, 

and for almost all (a,x) g [— r, 0] x [0, a], 

7n(M k M k N k (f 9 (f,g)))(<r,x)>0. 

Then the result is the same. If we assume that g(x) > 0 for almost all x G [xo,xi] with some 
0 < xo < x\ <a. Then for almost all x G [xo, a], 

7t 2 [7t 2 (N k (f 9 (f 9 g)))](x)>0. 
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This further implies that for almost all x e [xq, a], 

7T 1 [7t 2 (M k N k (f,(f,g)))](x)>0. 

Then for almost all x e [0, a], 

TtdMMiMMtf, (f,g)))](x) > 0, 7t 2 [7t 2 (M k M k N k (f, (f,g)))](x) > 0, 

and for almost all (cr,x) e [— r, 0] x [0,a], 

7n(M k M k M k N k (f, (f,g)))(cr,x) > 0. 

Then the result is also the same. If we assume that v = 0 or v = 1, the arguments are the similar 
by the assumption (H.4). This completes the proof. ■ 

By the positivity of the semigroup (T(7))?>o> s(A) > — oo(see Theorem VI. 1.10 of Engel and 
Nagel [8]), Corollary 3.2 and Lemma 3.3, we have the following result: 

Lemma 3.5 s(A) e a (A) and a; ess (A) < s(A). 

By Lemmas 3.1-3.5 and Corollary V.3.3 of Engel and Nagel [8], we conclude that there exist 
a positive rank one projection operator n in E and constants s > 0 and M > 1 such that 

||e-* (A) 'T(7) - n|| <Me~ £t fort> 0, (21) 

where || • || denotes the operator norm in E. This completes the proof of Theorem 1.2. 



4. Relation with the one-phase model 

In this section we give the proof of Theorem 1.3 and give a comparison with the sum of the 
densities of individuals in the 'reproductive' stage and the 'nonreproductive' stage N(t,x) and 
the solution N(t,x) of the problem (6). The idea is inspired by the generalized relative entropy 
method [12,13]. In [3], we use the similar method to give a comparison between a two-phase cell 
division model with the classical one-phase model. 

Proof of Theorem 1.3 By Lemma 3.5, we know that A = s(A) is the dominant eigenvalue of the 
eigenvalue problem 

a 

Xw(cr,x) w(cr,x) =0, — r < o < 0, 0 < x < a, 



da 
dx * 



Xu(x) + -j-(y(x)u(x)) = — /jLi(x)u(x) — p\(x)u(x) + p2(x)v(x) 



na /»0 



+ v I I P(cr,x,y)w(a,y) da dy, 0 < x < a, 



d (22) 
Xv(x) + —(y(x)v(x)) = -fi 2 (x)v(x)pi(x)u(x) - p 2 (x)v(x) 
dx 

pa pO 

+ (1 — v) / / fi(a,x, y)w(a, y) da dy, 0 < x < a, 

w(0,x) = u(x), 0 < x < a, 
fi(0) = 0, v(0) = 0 
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and the corresponding eigenvector (w(a,x), (S, v)) is strictly positive, i.e. w(a,x) > 0 for all 
(cr,x) e [— r, 0] x (0, a), u(x) > 0 and v(x) > 0 for all 0 < x < a. Letting 

U(x) = u(x) + v(jc), 0(x) = U{X) 



U(x) 

fi(x) = 0(x)fii(x) + (1 - 0(x))ii 2 (x), Pi(cr,x,y) = 0(x)P(cr,x,y), 

we have that 

d ~ 

XU(cr,x) U(cr,x) = 0, — r < o < 0, 0 < x < a, 

d(y(x)U(x)) « r r° 

XU(x) -\ = — /jl(x)U(x) + / / #i (a, x,y)U (a, y) da ay, 0 < x < a, 

dx Jo J-r " " (ZJj 

£/(0,jc) = U(x), 0 <x <a, 
U(0) = 0. 

Hence, X = 5(A) is also an eigenvalue of the above eigenvalue problem, and (U(a,x), U(x)) is 
the corresponding eigenvector. Since U(a,x) > 0 for all (a, x) e [— r,0] x (0,a) and > 0 
for all 0 < x < a, it follows that X = s(A) is also the dominant eigenvalue of the above eigenvalue 
problem. 

Now let (p(t,x), n(t,x)) and N(t,x) be the solutions of the model (3) and (5) with initial data 
(p(a,x),no(x)) and N (a, x) = p(a,x) + «oW, respectively, and set 

N(t,x) = p(t,x) + n(t,x) for f > 0, 
N(cr,x) = p(cr,x) + no(x) for cr e [-T, 0]. 

We want to compare N(t,x) and N(t,x). From Equation (21) and Theorem 9.1 1 of Clement et ah 
[7], we have the following asymptotic expression: 

e~ Xt (p(t + a,x),(p(t,x),n(t,x))) = ci(w(a,x), (u(x), v(x))) + 0(e~ £t ) as t -> 00, (24) 

where c\ is a constant uniquely determined by the initial data (p, no); c\ provided (p(a,x), no 
(x)) > (0,0) (i.e. (p,n 0 ) > (0,0) and (p,n 0 ) # (0,0)). This implies that 

e~ Xt N(t, x) = ci U(x) + 0(Q~ £t ) as 00. (25) 

From Equation (21) of Bai and Xu [4] and Theorem 9.11 of Clement et al [7], we have the 
following asymptotic expression: 

e~ Xt (N(t + cr,x),N(t,x)) = c 2 (U(a,x), U(x)) + 0(t~ £lt ) as 00, (26) 

where C2 is a constant uniquely determined by the initial data N\ c 2 > 0 provided N > 0 (i.e. 
TV > 0 and TV 7^ 0). This also implies that 

e~ Xt N(t,x) = c 2 U(x) + 0(Q~ £2t ) as f 00. (27) 

From Equations (25) and (27), we obtain 

e~ Xt N(t,x) - e~ Xt N(t,x) = c 3 U(x) + 0(z~ £t ) as 00, (28) 

where C3 = c 2 — ci and e = minimi, £2}. In what follows we prove that, generally speaking, 
C3 / 0. Let (0, (cp, VO) be the eigenvector of the conjugate problem of Equation (22) (see Theorem 
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8.17 of Clement et al. [7]), i.e. 
A0(c 

OCT 

+ (1 — v) / fi(cr,y,x)\lf(y) dy, — x < o < 0, 0 < x < a, 



d r 

(a,x) + — 0(a,x) = v / fi(cr,y,x)(p(y)dy 
do Jo 

> ["Pier, 
Jo 



dcp(x) 

X(p(x) — y(x) = — ii\(x)(p(x) — p\(x)(p(x) + p\(x)^r(x) + 0(O,x), 0 < x < a, 

dx 

d\J/(x) 

Xi/f(x) — y(x) = — ii\(x)^r (x) + p2(x)^r(x) — p2(x)cp(x), 0 < x < a, 

dx 

0(— r,x) =0, 0 < x < a, 
<p(«) = 0, 
^(a) =0. 

(29) 

We normalize (0, (<p, i/>0) such that 
r0 



w(a,x)0(o-,x) dadx + / (S 

-r JO 



w(cr,x)(f)(cr,x) dadx + / («(x)^(x) + v(x)^ (x)) = 1. (30) 



By Theorem 8.17 of Clement et al. [7], 0(o%x) > 0 for all (a,x) e [— r,0] x (0, a), cp(x) > 0 
and ^(x) > 0 for all 0 < x < a, due to a similar reason as that for w, u and v. Next, let (O, ^) be 
the eigenvector of the conjugate problem of Equation (23) (see Theorem 8.17 of Clement et al. 
[7]), i.e. 

d r 

AO(a,x)H <I>(cr,x) = / ft(cr, j,x)vl/(y) dv, — r < a < 0, 0 < x < a, 

dor J 0 

(M(jc) 

A*I>(x) - y(x) = -fi(x)W(x) + O(0,x), 0 <x <a, (31) 

dx 

0(-r,x) =0, 0 < x < a, 
vl/ (a) = 0. 

We normalize (O, ^) such that 

U(cr,x)<$>(a,x)dcrdx+ / U(x)V(x) dx = 1. 

-r JO 

By Theorem 8.17 of Clement efaZ. [7],0(a,x) > Ofor all (a,x) e [-r,0] x (0, a) and *J> (x) > 0 
for all 0 < x < a. From Equations (3) and (29), we easily obtain 



<S (/ / /7 ^ + or ' Jc ^^ a '^ e "^dadx + y 0(/,jc)^(jc) + /i(f,Jc)^(Jc))e _A,r djc^ = 0. 
Hence 

n0 /»a 
+ a,x)0(o-,x)e _A? do'ck + / (p(£,x)<p(x) + jc)^t(jc)) e" A? dx 
-r JO 

n0 
/?(<7,x)0(a,x) dadx + / (poW^W + ^o WV^ W) dx 
-r Jo 
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for all t > 0. Letting t -> oo and using (24), we get 



na nO na 

c\ / w(a,x)(/)(a,x) da dx + / + v(x)ty (x)) dx 

J0 J-t Jo 

na nO na 

/ / p(a,x)4>(a,x) da dx + / (po(x)(p(x) + no(x)ty(x)) dx. 
Jo J-t Jo 



From Equation (30), we obtain 

ra nO 



na nO na 

d = I / p(a,x)4>(a,x) da dx + / (po(x)(p(x) + no(x)\ls(x)) dx. 

Jo J-t Jo 



By a similar argument we have 

r a nO 



na nO na 

c 2 = / N(a,x)®(a,x)dadx+ / (p 0 (x)^ (x) + n 0 (x)^ (x)) dx 

Jo J-t Jo 

na nO na nO 

= p(a,x)$>(a,x)dadx+ / (po(x)^ (x) + n 0 (x)(^> (x) + / 0(a, jc) da)) dx, 

JO J-T JO J-T 



so that generally speaking we have c\ ^ c 2 or c 3 / 0. This proves Theorem 1.3. ■ 
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